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We present a condition for automatic O(a) improvement in twisted mass lattice QCD, using sym- 
metries of the Symanzik effective theory. If the continuum part of the Symanzik effective theory 
is invariant under a particular transformation, named T\ in this report, scaling violations of all 
quantities invariant under T\ transformation are even in the lattice spacing a. On the other hand, 
quantities non-invariant under T\ vanish in the continuum limit with odd powers in a. We prove 
this statement even for the massive case without using the equation of motion. We also consider 
a few different criteria for the T\ invariant condition in lattice theories and discuss ambiguities of 
the lattice condition for 0(a) improvement. 
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1. Introduction 

It becomes more and more apparent that twisted mass Lattice QCD (tmLQCD) [[T], is a 
promising formulation to approach the chiral limit of QCD, despite the fact that the flavor symme- 
try is explicitly broken. A twisted mass protects the Wilson-Dirac operator against small eigenval- 
ues and therefore solves the problem of exceptional configurations |J], thus making numerical 
simulations with small quark masses feasible [g]. This numerical advantage of tmLQCD is supple- 
mented by the property of automatic 0(a) improvement rfq]. For a recent review of these and some 
more results in twisted mass LQCD see Ref. [^]. 

The so-called "maximal twist" condition, required for the proof of automatic 0(a) improve- 
ment [g], however, causes some confusions. For example, it has been suggested [^] that maximal 
twist can be achieved by setting the bare untwisted mass to the critical quark mass of the Wilson 
fermion where the pion mass vanishes (we call this choice "the pion mass definition" in the fol- 
lowing). However, it has been pointed out ^ that this choice does not lead to automatic 0(a) 
improvement. Indeed, terms linear in a and with fractional powers of a are predicted by Wilson 
Chiral Perturbation Theory (WChPT) for very small twisted quark masses. On the other hand, 
automatic 0(a) improvement is expected to hold if the critical mass is defined through the partially 
conserved axial vector Ward identity quark mass (PCAC mass definition). 

In this report we present an explicit condition based on symmetries of the Symanzik theory, 
required for automatic 0(a) improvement, and prove that scaling violations of all quantities which 
have non-zero values in the continuum limit are even in a (0(a) improvement). The detailed proof 
of this statement has already been published in Ref. [ph. Therefore, in this report we try to avoid 
unnecessary rigorousness in our proof and stress the mechanism which leads to 0(a) improvement 
in twisted mass QCD. 

2. Condition and proof for automatic 0(a) improvement 

We first give our statement, which will be proven in this report. Twisted mass QCD with a 
certain condition leads to automatic 0(a) improvement, which means that operators as well as the 
action are automatically 0(a) improved without any improvement coefficients. This holds true 
even in the massive case and without the explicit use of the equations of motion. More explicitly, 
all scaling violations of non-zero physical quantities are even in a, while quantities which vanish 
in the continuum limit have only odd powers in a. 

2.1 Main idea of the proof 

The twisted mass lattice QCD action for the 2-flavor theory is given by S im = Sg + Sf, where 
So is the gauge action and 

Sf = J>z.(x)i [Y^V+ + V-) VL -arV+V- VL ] (x) +£ y L (x)M e i0 ^ 3 y L (x) (2.1) 

x,n x 

is the 2-flavor Wilson fermion action with a twisted mass term, where Mo and do denote the bare 
mass and bare twist angle. It is also customary to write Moe' 9 " 75 * = mo + //io/s^ 3 , using the bare 
untwisted mass m and the bare twisted mass jio- 
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This action is invariant under the following global transformations: (1) U(1)<8>U(1) vector 
symmetry, y L — > e l ( a o+ a 3 t ) \j/ L , y L — > )// L e _ '( ao+0!3 ' r ) . This transformation is part of the U(2) flavor 
symmetry of the untwisted theory. (2) Extended parity symmetry Pp' 2 : — ► T l ' 2 j4yi(Px), 
Wl(x) — > Yl(Px)Y4^ 1,2 where P is the parity transformation. Alternatively, one can also augment P 
with a sign change of the twisted mass term jUo, P = P X [jUo —>■ — jUo], which is also a symmetry of 
the action. (3) Standard charge conjugation symmetry, as in the untwisted theory. 

The lattice theory can be described by an effective continuum theory (the Symanzik theory), 
whose effective action is restricted by locality and the symmetries of the underlying lattice theory. 

Taking into account the symmetries listed above one finds S e g = So + aS\ +a 2 S2 H , where the 

first two terms are given as 

So = S „ + J d 4 xy(x) [y^ +M R e ie ^ 3 ] y(x), (2.2) 
Si=C x j d A xy(x)o^ v F^ v (x)y(x). (2.3) 



5o, gaU ge denotes the standard continuum gauge field action in terms of the gauge field tensor F^ v . 
The second term in So is the continuum twisted mass fermion action. It is worth mentioning 
that there is no "twisted" Pauli term yysT 3 o^vF^y present in Si, since such a term violates the 
symmetry P. 

In addition to the effective action we have to specify the direction of the chiral condensate, 
since chiral symmetry is spontaneously broken. From the fact that the direction of the chiral con- 
densate is completely controlled by the direction of the symmetry breaking external field (i.e. the 
quark mass) in the continuum theory, we can take 



v(M R ) 



J i 



(2.4) 

pa 



where limM^ohrny^oV^s) ^ 0. Here the vacuum expectation value (VEV) is defined with 



respect to the continuum action So- To say it differently, the VEV ( |2.4[ ) defines the twist angle 8 in 
the Symanzik theory. 

We now want to argue that the choice 6 = n/2 (or — Tt/2) corresponds to "maximal twist". In 
terms of the mass parameters this is equivalent to Mr = ^Lr and rtiR = 0. In this case the action and 
the VEVs become 

So = So,gauge + j d A x y(x) [y^D^ + IMrYsT 3 } y(x), (2.5) 
(YY)so = 0, (i//7y5T 3 y/)s = v(M R ) . It is easy to check that So, the continuum part of the ef- 

12 _ ■ 12 

fective action is invariant under y — > e lwr5Z y — > ye lw ^ iX , and therefore also under the Z2 
subgroup T\ of this continuous transformation, defined by T\\y = iy^x^yf, Tiy = xjfiysl : . Since 
Tj 2 = 1 in the space of fermion number conserving operators, which contain equal numbers of 
Y and y, the eigenvalues of T\ are 1 (7\-even) or —1 (7\-odd). The crucial observation is that 
the VEVs (yy) and {yijsT'y) are also invariant under this transformation. The T\ symmetry 
is not spontaneously broken, hence it is an exact symmetry of the continuum theory. The 0(a) 
term aS\ = aC\ J d 4 x y(x)a jlv F^ v (x)y(x), on the other hand, is odd under 7\. Therefore, non- 
vanishing physical observables, which must be even under T\, can not have an 0(a) contribution, 
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since the 0(a) term is odd under T\ and therefore must vanish identically. This is automatic 0(a) 
improvement at "maximal twist". Note that non-invariant, i.e. Ti-odd quantities, which vanish in 
the continuum limit, can have 0(a) contributions. 

The above argument gives just the main idea of our proof for automatic 0(a) improvement, 
and we will give a detailed proof in the next subsection. However, one of the most important points 
of our analysis is that the condition for automatic 0(a) improvement is the invariance of theory 
under T\ transformation, or more generally its continuous version, which corresponds to a part of 
the exact vector symmetry in continuum QCD at "maximal twist". One might even say that the 
T\ invariance is more fundamental for automatic 0(a) improvement than the notion of "maximal 
twist", which is one of the consequences of 7\ invariance. 

2.2 General proof 

Let us consider an arbitrary multi-local lattice operator 0[ P f({x}), where {x} represents 
x\ ,X2, ■ ■ ■ ,x n , d is the canonical dimension of the operator, t = 0, 1 and p = 0, 1 denote the transfor- 
mation properties under T\ and P: 

Ti : 0\lf({x}) - (-\yo\lf({x}), P: 0\lf({x,t}) - (-l)P \ P f {{-x,t}). (2.6) 

Here we do not include the dimension coming from powers of the quark mass in the canonical 
dimension d of the operators. 

The lattice operator 0[ P { corresponds to a sum of continuum operators O t,,p "' n (n: non-negative 
integer) in the Symanzik theory as 

oo 

<=Z^I4/*" (2-7) 

n=d t„,p„ 

where n is the canonical dimension of the continuum operator O tnPn,n which consists of \jf, \ff, 
and only without any mass parameters, and 

7i : 0'" p " n ({x}) -> (-l)'"<9'" p -"({x}), P : 0'" p "' n ({x,t}) -» (-l) p "0'" p - n ({-x,t}),(2.S) 

with t n ,p„ = 0, 1. To have a total dimension d in the expansion in Eq. (fO|), the coefficients c?t n 
must be dimensionless. Here, to make an argument simpler, we consider the lattice operator, whose 
power divergences can be subtracted without spoiling our proof [QJ. 

The following selection rules among these operators are crucial for our proof of automatic 
0(a) improvement: 

t + p + d = t„+p„ + n mod(2), p + = p n + (#^o)„ mod(2), (2.9) 
where #;Uo and (#jlio) n denote the numbers of /^o's in 0[ P f and c^ %n , respectively. The second 



equality can be easily proven by the invariance of the lattice action (gjj) under the P = P x [/ii 



•H) 



-jUo] transformation. To prove the first equality, we introduce the following transformation: 

(2.10) 



U:< 



(A^x) -» -A^(-x) 
Xj/(x) -» (e i7C ^f /2 yr(-x) 



e 



,3/2 
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which is a modified version of the transformation ^ introduced in Ref. [Qj. Since it is easy to 
show that the lattice action (2T) is invariant under T\ x <^J, in addition to the invariance under P^, 
the lattice action is invariant under T\X S! l d xP F . On the other hand, we can easily see that f^j x P^ 
counts the canonical dimension times the parity of the operator as 

9\ *Pf : Otf({x,t}) - (-l^'O^tf* ,-*}), (2.11) 
®\ xP F : 0^»>" ({*,*}) (-1)" +P "0 to '"({x,-*}). (2.12) 

Therefore, the invariance of the action under T\ x Q)\ x P F implies the first selection rule. 

Let us show how these selection rules are used to determine the structure of operators in the 
Symanzik theory. We first consider the Symanzik expansion of the lattice action S lm : 



Stm = £ [a 2 "(5 00 ' 2 "+ i Uoa5 11 >)+a 2 "- 1 ( J Uo^ 01 - 2 "- 1 +5 10 >- 1 )] (2.13) 
»=o 

DO 

= ^ + mSl l + ^[a 2n ^ n + a 2n -'s\ n _ l ] , (2.14) 

n=\ 

where S tnPn,n denotes the action in the Symanzik theory whose canonical dimension and transfor- 
mation properties under T\ and P are (d„,t„,n). Here we pull out the \La factor from terms which 
have odd powers in \x^a. Therefore remaining factors always have even powers in jj$a. To derive 
the first equality we use the selection rules such that = n + t„ + /?„mod(2) = p n + (#jUo) n mod(2) 
since the lattice action satisfies d + 1 + p = and p + #/io = 0. In the second equality we define 

S° 2n = ^ 2n + n S 01 - 2 "-\ mSl^S 10 '- 1 /^ =5 10 ' 2 "- 1 +Ato5 11 ' 2 "" 2 («> 1), (2.15) 

and the superscript or 1 denotes the transformation property under T\ as evident from the above 
definition. Similarly we have 



0£m = °° + I V n °^2n + +a 2n - X 0\ +2n _ x ] , (2.16) 

n=l 

OO 

ObM = 0\ + £ [« 2 "Oi +2 „ + +a 2 "- 1 0° d+2n _ i ] , (2.17) 

for the multi-local operator with the canonical dimension d, where again the superscript or 1 
denotes the transformation property under T\ . 

Now we can specify the condition for automatic 0(a) improvement: It is stated that the con- 
tinuum part of the action d2. 14 ) is invariant under T\ . This condition leads to m = 0, so that the 
continuum part of the action is given solely by Sq. We will consider the scaling behaviour of the 
vacuum expectation value of an arbitrary multi-local operator, (0[ at d ({x})). For this purpose we 
define 

{OO I OQ 

£ [a 2n ^ 2n + a 2n - l S\ n _ l ] \ = e$ £ a "S {n \ (2.18) 
n=\ ) n=0 

where we define a"S^ to be the sum of the a" terms in eq. ( |2.18[ ). For example, the first few terms 
are given as tf® = 1, = S\ , and = S° + {S\) 2 /2\ . Under the T\ transformation, they 
behave as T\ : — > (— l) n 5^ . By expanding both action and operator, we have 

oo oo 

(OU{*})> = E «"<<*(«)>*.= I «" + '(^(W)5 (/) ) 5 o (2.19) 

n=0 n=l=0 
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where t n = n + t mod (2). In the second line the T\ invariance tells us that (0^({x})S^) s o = 
unless t n + 1 = t + n + l = mod (2). Therefore we have 

oo 2s+t 

(OlAM)) = I> 2m I (0:;'({x})5( 2 ^-»)) 5 o, (2.20) 

s=0 n=0 

from which we derive 

«t(M)> = (O°({x})) s o + O(a 2 ) + O(a*) + -- (2.21) 
«({*})> = 0(a) + 0(a 3 ) + 0(a 5 ) + ■■■. (2.22) 

This proves our statement that the scaling violation of all Ty invariant operators, which have non- 
zero VEV in the continuum limit, are even in a, while that of T\ non-invariant operators, whose 
VEV vanish in the continuum limit, are odd in a. This is true for non-zero jUo and does not require 
the use of the equation of motion. 

2.3 Condition for 0(a) improvement in the lattice theory 

In the Symanzik theory, the condition for 0(a) improvement is uniquely defined by the condi- 
tion that an arbitrary T\ non-invariant operator O t=ip,d has a vanishing expectation value. Provided 
this condition is fulfilled, the expectation values of all 7\ non-invariant operators vanish. Hence 
the particular choice for lp ' d is irrelevant, and in that sense the condition is unique. In the lattice 
theory, however, the condition defined by (0^ t ,d ) = depends on the choice of the operator o\^ d , 
and is therefore not unique. In terms of the Symanzik theory, for this condition to be satisfied, an 
equation 

aFo(a 2 ,ma, IM)) + mF\(a 2 ,ma, /lo) = 0, 

where Foj are some functions of a 2 , ma and /lo, must be fulfilled by tuning the untwisted mass 
m. Since this equation is invariant under (m,a) — > (—m,—a), the solution has the form that 
m = af(a 2 ,Ho) under the assumption that the solution to the equation is unique. If one takes a dif- 
ferent lattice operator to define the T\ invariant condition, the solution is given by m' = af'(a 2 ,pL). 
Therefore the difference between two definitions is 0(a): m — m' = a(f — f). Note that a solution 
m in general depends on }Xo, inherited from the /i dependence of Fq \. 

Let us consider some examples for the condition in the lattice theory. A simple one is given by 
((Y¥)iat) = ■ Unfortunately, this definition is not very useful in practice, since the subtraction of 
power divergences necessary for {yy) prevents a reliable determination of this VEV in the lattice 
theory. Instead one may take 0\ at (x,y) = A a ^(x)P a (y) or Oi at (x,y) = d iX A a ^ l (x)P a (y) {a = 1,2), as 



was done in Refs. UlQ] |11|, |12J : 

(A£(x)P*(y)) = or (d^(x)P"(y)) =0, (2.23) 
where A a ^ and P" denote the axial vector current and pseudo scalar density, respectively. Yet another 



choice is [13] 



(Al(x)P\y))=0. 

Depending on the choice for the axial vector current, either the local or the conserved one, the 
conditions ( |2.23| ) lead to a different definition for maximal twist. However, the difference will be 
again of 0(a). 
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We close this section with a final comment. Any condition for 0(a) improvement in the lattice 
theory determines a value for the bare untwisted mass niQ as a function of the bare twisted mass 
Ho. It has been suggested to tune the untwisted mass to its critical value mo = m cr where the pion 
mass vanishes in the untwisted theory. However, this condition is not related to T\ invariance. 
For example, contributions from excited states violate eq. ( |2.23[ ) even at m K = 0. Consequently, 
the pion mass definition does not correspond to automatic 0(a) improvement according to the 7\ 
invariance condition. 

3. Conclusion 

In this paper we gave a comprehensive proof for automatic 0(a) improvement in twisted mass 
lattice QCD. The most important observation is that a precise definition for 0(a) improvement is 
described by the symmetry in the continuum theory. If the continuum part of the Symanzik theory is 
invariant under T\ transformation, scaling violations for all quantities are shown to be even powers 
in a, as long as they are invariant under the T\ transformation. Non-invariant quantities, on the 
other hand, vanish as odd powers in a. 
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